Astronomy & Astrophysics manuscript no. tf 'lang 


April 8, 2008 


(DOI: will be inserted by hand later) 





IMAGES*-III: The evolution of the Near-Infrared Tully-Fisher 

relation over the last 6 Gyr 

M. Puech 1 - 2 , H. Flores 2 , F. Hammer 2 , Y. Yang 2 , B. Neichel 2 , M. Lehnert 2 , L. Chemin 2 , N. Nesvadba 2 , B. Epinat 5 , P. 

Amram 5 , C. Balkowski 2 , C. Cesarsky 1 , H. Dannerbauer 6 , S. di Serego Alighieri 7 , 1. Fuentes-Carrera 2 , B. 
Guiderdoni 8 , A. Kembhavi 3 , Y C. Liang 9 , G. Ostlin 10 , L. Pozzetti 4 , C D. Ravikumar 1 1 , A. Rawat 2 - 3 , D. Vergani 12 , J. 

Vernet 1 , and H. Wozniak 8 

1 ESO, Karl-Schwarzschild-Strasse 2, D-85748 Garching bei Miinchen, Germany 

2 GEPI, Observatoire de Paris, CNRS, University Paris Diderot; 5 Place Jules Janssen, 92190 Meudon, France 

3 Inter-University Centre for Astronomy and Astrophysics, Post Bag 4, Ganeshkhind, Pune 41 1007, India 

4 INAF - Osservatorio Astronomico di Bologna, via Ranzani 1, 40127 Bologna, Italy 

5 Laboratoire d'Astrophysique de Marseille, Observatoire Astronomique de Marseille-Provence, 2 Place Le Verrier, 13248 
Marseille, France 

6 MPIA, Konigstuhl 17, D-69117 Heidelberg, Germany 

7 INAF, Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5, 1-50125, Florence, Italy 

8 Centre de Recherche Astronomique de Lyon, 9 Avenue Charles Andre, 69561 Saint-Genis-Laval Cedex, France 

9 National Astronomical Observatories, Chinese Academy of Sciences, 20A Datun Road, Chaoyang District, Beijing 100012, 
PR China 

10 Stockholm Observatory, AlbaNova University Center, Stockholms Center for Physics, Astronomy and Biotechnology, 
Roslagstullsbacken 21, 10691 Stockholm, Sweden 

11 Department of Physics, University of Calicut, Kerala 673635, India 

12 IASF-INAF - via Bassini 15, 1-20133, Milano, Italy 



Received ; accepted . 



Abstract. Using the multi-integral field spectrograph GIRAFFE at VLT, we have derived the K-band Tully-Fisher relation 
(TFR) at z ~ 0.6 for a representative sample of 65 galaxies with emission lines (W (OII) > 15A). We confirm that the scatter 
in the z ~ 0.6 TFR is caused by galaxies with anomalous kinematics, and find a positive and strong correlation between the 
complexity of the kinematics and the scatter that they contribute to the TFR. Considering only relaxed-rotating disks, the scatter, 
and possibly also the slope, of the TFR, do not appear to evolve with redshift. We detect an evolution of the K-band TFR zero 
point between z ~ 0.6 and z = 0, which, if interpreted as an evolution of the K-band luminosity of rotating d isks, would imply 
that a brightening of 0.66±0. 14 mag occurs between z ~ 0.6 and z = 0. Any disagreement with the results of Flores et al. I fcoodh 
are attributed to both an improvement of the local TFR and the more detailed accurate measurement of the rotation velocities in 
the distant sample. Most of the uncertainty can be explained by the relatively coarse spatial-resolution of the kinematical data. 
Because most rotating disks at z ~ 0.6 are unlikely to experience further merging events, one may assume that their rotational 
velocity, which is taken as a proxy of the total mass, does not evolve dramatically. If true, our result implies that rotating disks 
observed at z ~ 0.6 are rapidly transforming their gas into stars, to be able to double their stellar masses and be observed on the 
TFR at z = 0. The rotating disks observed are indeed emission-line galaxies that are either starbursts or LIRGs, which implies 
that they are forming stars at a high rate. Thus, a significant fraction of the rotating disks are forming the bulk of their stars 
within 6 to 8 Gyr, in good agreement with former studies of the evolution of the mass-metallicity relationship. 

Key words. Galaxies: evolution; Galaxies: kinematics and dynamics; Galaxies: high-redshifts; galaxies: general; galaxies: 
interactions; galaxies: spiral. 



1. Introduction 

Since the firs t rotation curves were measured at intermedi- 
ate redshifts dVogt et alj 1 1993b . many studies have been de- 
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voted to the evolut ion of the Tully-Fisher Relation (TFR, 



Tullv & Fisher 1 1977b . given it s prominent role i n constraining 



galaxy-formation models (e.g jDutton et al ] l2007l) . Early work 
usin g B-band imaging foun d only modest luminosity evolu- 
tion ( Vog t et al . 1996, 1997), but subsequent studies using data 
of the same band have suggested a more significant evolution 
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(Simard & Pritchet 1998; iBamford et alj 120061; iWeiner et al 
20061; IChiu et al.ll2007l)~ In almost all of these studies, the B- 



band TFR shows a large dispersion in comparison with the lo- 
cal re lation, especially a t the low-luminos i ty (or velocity) end 
(e.g., lBohmet al.ll2004 . iBohm & Zieglerl d2007l) showed that 
this effect could be attributed to an incompleteness in magni- 
tude provided that the scatter decreased by a factor of at least 
three between z ~ 0.5 and z — 0, with no evolution of slope or 
zero point. Alternatively, they proposed a possible luminosity- 
dependent evolution, in which distant low-luminosity galaxies 
would have lower mass-to-light ratios [M/L] than their local 
counterparts, while higher-luminosity galaxies would not un- 
dergo strong M/L evolution. Therefore, the situation for B- 
band data remains unclear, particularly concerning the large 
scatter in measurements found at high redshift. 

To minimize observational biases and enable easier com- 
parisons with models, a progressive interest has been de- 
voted to the K-band and stellar-mass TFRs (hereafter, smTFR). 
Long-slit spectroscopy of distant galaxies revealed an smTFR 
of significant scatter, with detectabl e evolution of neithe r 



zero point nor slope, up to z ~ 1 dConselice et al.l 120051) 



This much larger dispersion appears to be linked to "kine- 
mati cally anomalous galaxies", as inferred from local stud- 
ies ( Kannappan & Bartonl 12004 ). or, from a morphological 
point of view, with disturbed, compact, or merging galaxies 
dKassin et all2007t lAfkinson et alj[2007b . IWeiner etalj d2006h 
and iKassin et al.l d2007l) defined a new tracer of the galaxy 
dark-halo potential called S , which corrects the rotation veloc- 
ity for disordered, non-circular motions. Once expressed us- 
ing this new kinematical estimator, the distant smTFR shows 
significantly -reduced scatter, with no detectable evolution in ei- 
ther zero point or slope. This suggests an important role of non- 
ordered motions in increasing the scatter of the distant TFR. 

In this respect, 3D spectroscopy provides a unique way 
to distinguish relaxed rotating disks from other kinematically- 
disturbed galaxies. Kinematically-selected relaxed (or pure) ro- 
tating disks at z ~ 0.6 present a TFR tha t appears to evolv e 
neither in slope, zero point, nor dispersion ( Flores et alj|2006l) : 
kinematically-anomalous galaxies, which range from perturbed 
rotators where a rotation can be detected to galaxies with com- 
plex kinematics but no noticeable rotation, appear to be respon- 
sible for the increased scatter of the distant TFR. Therefore, 
3D spectroscopy allows us to establish a direct connection be- 
tween non-ordered motions and kinematical and morphologi- 
cal anomalies In this paper, we double the sample previously 
available (see lYang et al.l l2007[ hereafte r Pape r I) to test ro- 
bustly the results obtained in lFlores et al. I d2006l) . 



This paper is organized as follows: Sect. 2 introduces the 
data used in this study; Sect. 3 details the methodology used 
for the analysis; Sect. 4 presents the K-band TFR obtained 
from GIRAFFE data, and Sect. 5 draw the conclusions from 
this work. In two appendixes, we derive the K-band luminos- 
ity density and the stellar-mass TFR in the GIRAFFE sample. 
Throughout, we adopt Ho = 70 km/s/Mpc, Q.m = 0.3, and 
Qa = 0.7, and the AB magnitude system. 



2. Data 

2.1. Distant sample 
2.1 .1 . Kinematics 

We used the multi-object integral field spectrograph FLAMES - 
GIRAFFE at VLT, to obtain the [Oil] spatially-resolved kine- 
matics of a sample of 68 galaxies at 0.4 < z < 0.75. The 
median redshift of the sample was z me d ~ 0.61, which cor- 
responds to a look-back time ~ 5.8 Gyr, while the 25 and 
75 percentiles of the redshift distribution are Z25 ~ 0.49 and 
Z75 ~ 0.67. This sample represents well the luminosity func- 
tion of z ~ 0.6 galaxies with stellar masses [M sle u ar ] in the 
1 .5 - 15x10 1() Mq range, and is unaffected by field-to-field vari- 
ations within Poisson statistics, as targets were observed in four 
dif ferent fields ( s ee P aper I). 

d2006l) and Paper I used GIRAFFE velocity 
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fields [VF] and velocity dispersion maps [cr-maps], to de- 
fine three kinematical classes, namely the rotating disks [RD], 
which correspond to relaxed rotators, the perturbed rotators 
[PR], corresponding to rotating disks showing some weak dis- 
turbances, and kinematically -complex [CK] galaxies, corre- 
sponding to dynamically non-relaxe d galaxies, probab l y asso - 
ciated with major mergers (see also IPuech et al. 2006 , 2007 ). 



This classification takes into account the residuals between the 
observed VF and <x-map and those predicted by a rotating-disk 
model (see Paper I), which mitigates the subjectivity of a fully 
visual classification. 

2.1.2. Morphology 

A detailed morphological analysis of a subsample of 52 galax- 
ies, whic h had multi-band HS T/ACS imaging data, was com- 
pleted bv lNeichel et al.l(l2008l) (hereafter, Paper II). They found 
a good agreement between both kinematical and morphological 
classifications. Only 16% of galaxies in this subsample are both 
kinematically-classified as RDs and morphologically-classified 
as spirals. These "rotating spiral disks" were selected to have 
a redder bulge than the disk and to be therefore quite similar 
to local spirals, apart their much bluer integrated colors and 
[Oil] equivalent widths (see Paper II for a detailed discussion). 
Furthermore, they showed no special trend in size nor in Bulge- 
to-Total luminosity ratio [B/T] compared to local galaxies, and 
we do not expect any bias in the distant TFR that could be 
due to morphologi cal variation o f the TFR along the Hubble 
Sequence (see, e.g.. lRussellll2004l) . In the following, we use the 
inclinations [i] and half light radii [Rhaif] derived in Paper II. 

The comparison between kinematics and morphology re- 
vealed two special cases of RDs. First, J033230.78-275455.0 
was detected in emission only on one half of the galaxy, which 
was caused by the superposition of a skyline on an [Oil] emis- 
sion line. This galaxy was classified as a RD, but the amplitude 
of its VF is affected by large uncertainty. Second, J033241.88- 
274853.9 was classified as a Tadpole galaxy in Paper II because 
of its highly-asymmetric shape. In the following (see Sect. 3.2), 
we assume that all RD flux distributions can be approximated 
by an exponential disk, which, obviously, does not apply to this 
galaxy. Therefore, in this paper, "RDs" refers to galaxies classi- 
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fied as rotating disks in Paper I excluding J033230. 78-275455.0 
and J033241. 88-274853.9, while "RD+" refers to all RDs in- 
c/z«#ngJ033230.78-275455.0and J033241. 88-274853.9. 



corrected for extinction), (10) Absolute K-band magnitude M C K 
corrected for extinction (±0.2 mag, see Sect. 2.1.3), and (11) 
Stellar mass in log (M, fe // ar /M ) (±0.3 dex, see Appendix A). 



2.1.3. Absolute K-band magnitudes 



Abso lute K-b and magnitudes M K w ere taken from lFlores et al 
(2006) and iRavikumar et al.l d2007l) They were derived using 
Bruzual & Chariot 2001 stellar-population models, assuming 
a CSP template with solar metallicity and an exponential star- 
formation history with r = 0.5 Gyr, which desc ribes the prop 



erties of most galaxies between z = 0.4 and 1 dHammer et al 



2001). For each galaxy, an optimal SED was found by match- 
ing the observed J-K color. This method was preferred to com- 
plete SED fitting because 50% of the selected galaxies, from 
the CFRS and the HDFS, did not have similar multi wavelength 
photometric data to galaxies in the remaining 50%, which were 
selected from the CDFS; SED fitting would, therefore, not have 
been able to measure an Mk of similar quality for all galaxies 
studied. We compared, for galaxies from the CDFS, the im- 
pact of this simple "color-matching" method relatively to a full 
SED fitting (using optical photometry only): we found an ex- 
cellent agreement between J-band absolute magnitudes with a 
1 -<x dispersion ~ 0.01 mag, and no systematic effect (Hammer 
& Pozzetti, private communication). We compared absolute J- 
band and not K-band magnitudes because whatever the method 
used, these are not extrapolations. A potential drawback of the 
"color-matching" method is that Mk is extrapolated, since the 
reddest photometric point used is the observed K-band, which 
falls roughly in the middle of the rest-frame J and H band at 
Z ~ 0.6. For galaxies in the CDFS, IRAC 3.6/im fluxes were 
publicly available, which allowed us to compare the Mk ob- 
tained using this method with those obtained taking into ac- 
count the IR 3.6/mi flux: we found a residual 1 - <x dispersion 
of 0.12 mag between both estimates, with no systematic effect. 
Using two independent photometric datasets for galaxies in the 
CDFS (i.e., EIS and ACS), we derived a random uncertainty of 
< 0.2 mag on Mk, with no noticeable systematic uncertainties. 
Therefore, we adopted a random uncertainty of 0.2 mag and 
neglected systematic uncertainties. 

Absolute magnitudes were correcte d for internal extin ction, 
using the mass-dependent method of [Tullv et al.l (11998). We 
applied an additional correction of 0.04 mag that accou nted 
for extinction in face-on galaxies (Tully et al. 1998, 1985). We 
note that within the sample of 68 objects, three do not have NIR 
photometry (see TableQ]). 

2.1 .4. Summary 

For the distant sample, all measurements used are listed in 
Table [1] The columns of the table are as follows: (1) IAU ID, 
(2) Internal ID, (3) Redshift z, (4) Inclination i (in deg., ±5 
deg, see Sect. 2.1.2), (5) Dynamical Class D.C., (6) Projected 
radial velocity AV b s sin (/) (in km/s, see Sect. 3.1), (7) Rotation 
velocity corrected for inclination and spatial-resolution effects 
Vfiat (in km/s, see Sect. 3.1), (8) Total uncertainty in V/i at (in 
km/s, see Sect. 3.4), (9) Absolute K-band magnitude M K (un- 



2.2. Local sample 

As a local reference , we adopt the K-band TFR derived by 
(2007) for a complete subsample of the SDSS 
2007b . which allows us to control systematic ef- 
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(Pizagnoetal 



fects that can occur when comparing local and distant samples. 

One important choice for studying the TFR is the kinemat- 
ical estimator used for the rotation velocity V rot . Studies of the 
local TFR have shown that using different estimators (e.g., the 
maximal rotation velocity V„ wx , the plateau rotation curve ve- 
locity Vfiat, the velocity measured at the radius containing 80% 
of the light Vgo, and the velocity V2.2 measured at 2.2 disk scale 



length) can lead to different results (see, e.g.. lVer heiien 2001 
Pizagno et al ]l2007l) . Vfi at nas been shown to be correlated with 



V rot , since it is less influenced by the bulge dynamics, which 
can produce a central "bump" w ith V mnr > Vfiat- This can then 



lead to a tighter TF relation (e.g., Verhe iienl2001 ), and a signif- 
icant improve ment in the linearity of the K-band relation at the 
high-mass end (Noord ermeer & Verheiienll2007 ) . For the local 
samp le, we adopted the Vgo measurements of (IPizagno et al 



120071) using arctan fits to the RCs. To limit uncertainties, we 
restricted the local sample to galaxies for which Vgo is a good 
pr oxy for Vfiat (i.eJPizagno et al ,l2007l flags 1 and 2), as shown 
bv lHammer et all d2007l) . 



Rotation velocities were corrected for inclination using es- 
timates derived from their morphological axis ratio. When 
well-resolved 2D kinematics is available, it is possible to de- 
rive the inclination directly from the fit of the VF. However, 
large differences (up to ~ 10 deg) can be found between such 
kinematically-derived inclinations and the ones inferre d from 
the morphological axis ratio (e.g., IChemin et al.ll2006l) . Since 
we do not have 2D kinematics for galaxies in the local sample, 
we use morphologically-derived inclinations exclusively. This 
provides us with homogeneous estimates for the local and dis- 
tant samples, since in the latter we also use such inclinations. 

Hammer et al. I d2007l) combined the IPizagno et al.l (2007) 



kinematic measurements with 2MASS Ks-photometry. We es- 
timated absolute magnitudes in the local sample following a 
similar method to that used for the distant sample, including 
corrections for extinctions (see Sect. 2.1.3). The 2MASS Ks 
filter has the advantage of being close to the ISAAC Ks fil- 
ter used in the distant sample. Both filters match well the K- 
band LC O filter designe d to est a blish the faint IR standard-star 
system of IPersson et al. ( 19981) dMason et alj|2007t ICarpenteri 
200ll) . We assumed that both filters are identical during the 



SED-fitting procedure, which did not introduce any system- 
atic effect. We note than 2MASS K-band magnitudes were con- 
verted int o AB magnitudes us ing Ks(AB) = Ks{Vega) + 1.85, 
following EantonetaD d2005l) . 



Since we are exploring the higher tail of the stellar-mass 
distribution and the TFR is highly sensitive to stellar mass 
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Table 1. Principle properties of the sample of 68 galaxies used in this study, ordered by increasing RA (see text). 



IAUID 


Internal ID 


z 


i 


D.C. 


AV„f, r sin (;') 


Vfla, 
Jlar 


<-» Jlar 


M K 


Mi. 

K 


l02 (Maellnr/Mn) 
ivy t> V 1 '-* stellar I 


J030225.28+001325.1 


CFRS030085 


0.6100 


71 


RD 


127 


170 


23 


-20.83 


-20.98 


10.21 


J030228.72+001333.9 


CFRS030046 


0.5120 


66 


RD 


148 


200 


29 


-21.68 


-21.82 


10.64 


J030232. 16+000639.1 


CFRS039003 


0.6189 


29 


RD 


110 


290 


79 


99.99 


99.99 


99.99 


J030238.74+000611.5 


CFRS031032 


0.6180 


37 


PR 


83 


200 


36 


-22.64 


-22.70 


11.00 


J030239.38+001327.1 


CFRS030523 


0.6508 


41 


CK 


63 


130 


38 


-21.55 


-21.61 


10.48 


J030240.45+001359.4 


CFRS030508 


0.4642 


38 


CK 


47 


100 


40 


-20.34 


-20.40 


10.00 


J030240.99+000655.4 


CFRS031016 


0.7054 


67 


CK 


15 


30 


25 


-21.24 


-21.25 


10.37 


J030242.19+001324.3 


CFRS030488 


0.6069 


41 


CK 


36 


70 


36 


-20.84 


-20.89 


10.20 


J030245.67+001027.9 


CFRS030645 


0.5275 


45 


CK 


89 


150 


38 


-21.34 


-21.42 


10.40 


J030246.94+001032.6 


CFRS030619 


0.4854 


27 


RD 


70 


210 


76 


-21.93 


-21.99 


10.63 


J030248.41+000916.5 


CFRS031353 


0.6340 


57 


RD 


174 


260 


41 


-22.42 


-22.54 


10.85 


J030249. 10+001002.1 


CFRS031349 


0.6155 


48 


PR 


175 


290 


46 


-22.92 


-23.01 


11.04 


J030252.03+001033.4 


CFRS031309 


0.6170 


71 


CK 


84 


150 


22 


-22.91 


-23.05 


11.03 


J033210.25-274819.5 


CDFS4301297 


0.6100 


69 


PR 


110 


150 


26 


-20.83 


-20.96 


10.29 


J033210.76-274234.6 


CDFS4402679 


0.4180 


26 


CK 


186 


550 


123 


-23.62 


-23.68 


11.44 


J033212.39-274353.6 


CDFS 3400803 


0.4230 


90 


RD 


105 


180 


22 


-21.55 


-21.85 


10.61 


J033213.06-274204.8 


CDFS3500001 


0.4220 


79 


CK 


90 


130 


22 


-20.56 


-20.73 


10.19 


J033214.97-275005.5 


CDFS3300063 


0.6680 


22 


PR 


62 


190 


90 


-22.45 


-22.50 


10.92 


J033217.62-274257.4 


CDFS3401109 


0.6470 


46 


CK 


131 


250 


43 


-21.15 


-21.24 


10.38 


J033219.32-274514.0 


CDFS3400329 


0.7250 


72 


CK 


116 


270 


29 


-21.16 


-21.36 


10.39 


J033219.61-274831.0 


CDFS3300651 


0.6710 


49 


PR 


82 


190 


33 


-20.90 


-20.98 


10.27 


J033219.68-275023.6 


CDFS3202670 


0.5610 


58 


RD 


161 


230 


33 


-22.30 


-21.42 


10.88 


J033220.48-275 143.9 


CDFS3202141 


0.6790 


63 


CK 


39 


70 


24 


-20.63 


-20.70 


10.18 


J033224.60-274428.1 


CDFS 34006 18 


0.5380 


65 


CK 


63 


90 


27 


-20.35 


-20.44 


10.08 


J033225. 26-2745 24.0 


CDFS 3400279 


0.6660 


60 


CK 


47 


80 


26 


-21.54 


-21.61 


10.56 


J033226.23-274222.8 


CDFS3401338 


0.6679 


76 


PR 


150 


200 


24 


-21.93 


-22.13 


10.72 


J033227.07-274404.7 


CDFS 3400743 


0.7390 


84 


CK 


71 


110 


21 


-20.95 


-21.15 


10.26 


J033228.48-274826.6 


CDFS3300684 


0.6697 


22 


CK 


23 


130 


66 


-21.68 


-21.73 


10.63 


J033230.43-275304.0 


CDFS2200433 


0.6460 


70 


CK 


219 


380 


29 


-21.68 


-21.90 


10.64 


J033230.57-274518.2 


CDFS2400243 


0.6810 


34 


CK 


41 


80 


46 


-22.86 


-22.91 


11.08 


J033230.78-275455.0 


CDFS2102060 


0.6870 


66 


RD+ 


90 


200 


43 


-21.82 


-21.96 


10.66 


J033231.58-274121.6 


CDFS2500322 


0.7047 


42 


RD 


81 


140 


41 


-20.51 


-20.57 


10.10 


J033232.96-274106.8 


CDFS2500425 


0.4690 


16 


PR 


31 


210 


117 


-20.13 


-20.17 


10.01 


J033233.90-274237.9 


CDFS2401349 


0.6190 


17 


PR 


46 


200 


107 


-21.82 


-21.96 


10.66 


J033234.04-275009.7 


CDFS2300055 


0.7030 


59 


RD 


104 


160 


29 


-20.50 


-20.60 


10.09 


J033234. 12-273953.5 


CDFS2500971 


0.6280 


32 


CK 


25 


110 


44 


99.99 


99.99 


99.99 


J033237.54-274838.9 


CDFS2300477 


0.6650 


31 


RD 


106 


230 


71 


-21.98 


-22.03 


10.70 


J033238.60-274631.4 


CDFS2301047 


0.6220 


60 


RD 


118 


210 


33 


-21.45 


-21.57 


10.53 


J033239.04-274132.4 


CDFS2500233 


0.7330 


43 


PR 


58 


130 


39 


-20.62 


-20.68 


10.14 


J033239.72-275 154.7 


CDFS2200829 


0.4160 


35 


CK 


16 


30 


36 


-20.94 


-20.98 


10.31 


J033240.04-274418.6 


CDFS2400536 


0.5223 


16 


CK 


108 


470 


214 


-21.95 


-22.00 


10.77 


J033241. 88-274853.9 


CDFS2300404 


0.6680 


67 


RD+ 


80 


120 


27 


-20.90 


-21.01 


10.27 


J033243.62-275232.6 


CDFS2200611 


0.6800 


71 


PR 


38 


60 


23 


-19.94 


-20.01 


9.86 


J033244.20-274733.5 


CDFS2300750 


0.7365 


39 


CK 


46 


170 


40 


-21.76 


-21.83 


10.62 


J033245. 11-274724.0 


CDFS2300800 


0.4360 


43 


RD 


95 


270 


44 


-22.03 


-22.11 


10.80 


J033248.28-275028.9 


CDFS 1202537 


0.4462 


81 


PR 


76 


110 


22 


-20.38 


-20.54 


10.09 


J033249.53-274630.0 


CDFS 1302369 


0.5230 


46 


PR 


74 


150 


39 


-21.00 


-21.07 
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(see, e.g., iMcGaughl 120051) . we restrict the K-ban d TFR to 
log (Vfiat) > 2.2 and find from lHammer et al.l d2007l) : 



Radius [kpc @ z=0,8l 
4 6 



M K (AB) = -6.54 ± 1.33 - (6.88 ± 0.57) x log(V /to ), 



(1) 



with a residual dispersion cr res - 0.38 mag. Since we used uni- 
form Mk error bars, we used direct fits to the TFR, i.e., with 
M K as a V// ar dependent variable (using the IDL MPFITFUN 
procedure of C. Markwardt, translated from the MINPACK-1 
packag43)- 



3. Methodology 

3.1. Derivation of the rotation velocity 

For each galaxy, we estimated the deprojected VF half- 
amplitude using AV obs = (V 1 ™ - V;™)/(2 sin©), where V 1 ™ 
and VfFfa are respectively the maximal and minimal values of 
the VF sampled by the IFU (see Table [T). 

Because of the influence of the relatively coarse spatial- 
resolution of the kinemat ic data, AV b s underestimates the tru e 
rotation velocity V ro t (see lFlores et al .l2006l:IPuech et alJ2 006). 
Rather th an applying a mean correction factor to the entire sam- 
ple as in lFlores et al] (120061) . we corrected ech galaxy individ- 
ually by modeling its data-cube. We used a method developed 
from that used bv lFlores et al. (2006) and in Paper I, to model 
the cr-map. Assuming that all galaxies are thin-rotating disks, 
we modeled their data-cube, from their observed VF and cr- 
map, in the following way. 

First, for each galaxy, we constructed a grid of rotation 
curves (RC) with V ro , spaced at 10 km/s intervals, which 
roughly corresponded to the typical uncertainty in AV b s (see 
Sect. 3.4). Because of a lack of spatial resolution, the precise 
shape of distant RCs remains largely unknown. Therefore, we 
chose a simple arctan defined as VoX(2/7r)xarctan(r/r,), which 
depends only on two parameters, i.e., the a symptotic veloc - 
ity Vo, and the "turnover" radius r, (see, e.g., ICourteaul[l997l) . 
This RC shape was used in Paper I to model the VF of each 
distant galaxy. In the RD+ subsample, such models provided 
good matches to the observed VFs, which demonstrated that 
such a RC shape is a reasonable choice. The same shape was 
adopted in the local sample (see Sect. 2.2), which provided us 
with homogeneous estimates of the rotation velocity in both the 
local and distant samples. We chose not to fit a wide range of 
values for r, because a visual inspection of VFs revealed that 
for almost all RDs, t he gradient of the RC fell inside a single 
GIRAFFE IFU pixel dFlores et alJliooj Paper I). We therefore 
explored only three fixed values for r, (see Fig. [TJ, which al- 
lowed us to investigate more extreme cases where the RC is 
relatively steep or, in contrast, relatively fiat. 

The arctan RC model cannot reproduce the central "bump" 
observed in some early-type galaxies (see, e.g., Fig. [T). Given 
the coarse spatial resolution of the GIRAFFE IFU (0.52 arc- 
sec/pix, i.e., ~ 3.5 kpc at z = 0.6) and the large size of the see- 
ing disk compared to the typical scalelengths of distant galax- 
ies, it might however be impossible to distinguish clearly be- 
tween Vmax or Vfiat in distant galaxies. It might be expected 




0.5 1 1.5 

Eadius [aresec] 

Fig. 1. Illustration of the AV^ correction method, for an early- 
type RC (upper panel, note that V max > Vfi at ) and a late-type 
RC (bottom panel). The blue thick curves represent the input 
RC of the test, the black curves the best models (see text), and 
the dash-lines the two alternative RCs with non-optimal r,. We 
note that typical r, values lead to velocity gradients on spatial 
scales that are much smaller than the typical seeing of ~ 0.8 
arcsec. 



that the asymptotic velocity Vo, corresponding to an arctan RC, 
is probably, in this case, an average value between V max and 
Vfiat- Because we rely on Vfiat as a kinematic measure of the 
rotation velocity in the local sample (see Sect. 2.2), we are nat- 
urally led to quantify all random and systematic uncertainties 
relative to Vfi a t- Finally, we note that even for late-type RCs, 
Vo ~ Vfiat only if the RC has a well-defined plateau and if this 
plateau can be observed in terms of spatial coverage or SNR. 

From each RC, a high-resolution data-cube was con- 
structed, assuming a simple Gaussian shape for the emission 
line. For simplicity, we did not include noise in the simula- 
tions, assuming that the global uncertainty can be derived by 
estimating separately the influence of the noise and other ef- 
fects (see Sect. 3.4). As a consequence, neglecting the doublet 
spectral nature of the [Oil] emission line should not impact 
severally the results. The velocity width was assumed to be the 
minimal value observed in the cr-map, and the emission-line 
flux was taken from an exponential flu x distribution, assuming 
Rd ~ Rhaif/1-6 jPersic & Saluccifl~99 1 ) and limiting the spatial 



extent to a radius 2R 



half 



3Rd ~ R op t, where Rd is the disk 



http://cow.physics.wisc.edu/~craigm/idl/mpfittut.html 



scale-length, and R opl is the optical radius. We chose to use 
Rhaif rather than direct Rd measurements because more than 
80% of galaxies in the sample were not classified as simple spi- 
rals but had more complex morphologies (see Paper II), which 
could lead to meaningless Rd values. In contrast, Rhaif provides 
us with a uniform and well-defined size-parameter, which can 
be safely converted into Rd for thin exponential disks: using the 
bulge/disk decomposition of the RD subsample done in Paper 
II, we compared Rd with R^if /1.6, and found a maximal dif- 
ference of ~ 0.1 arcsec, which is much s maller than the typi - 
cal seeing during the observations (see also lPuech et al . 2007). 
Thus, such an error in Rd has little influence on the model, and 
we chose not to explore this parameter. 
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Each high-resolution data-cube was convolved by a 0.8 arc- 
sec seeing, which corresponded to the median condition of the 
observations, and then rebinned to the GIRAFFE sampling, i.e., 
0.52 arcsec/pix. From these simulated GIRAFFE data-cubes, 
simulated VFs and simulated AV mol j e i were derived as for real 
GIRAFFE data. We checked the influence of changing the see- 
ing from 0.8 to 1.0 arcsec, on the AV mo d e f- using Monte-Carlo 
simulations of 100 GIRAFFE data-cubes (see next section), we 
found a good linear correlation between the AV m „dei obtained 
for a 1 arcsec seeing and that obtained for a 0.8 arcsec seeing, 
all other properties being equal. This fit indicated that for a 1 
arcsec seeing, AV mo d e i was reduced by ~ 0.05% compared to 
that obtained for a 0.8 arcsec seeing (see Fig. [2j. Therefore, 
the maximal uncertainty on AV mot j e i due to seeing variations is 
~ 12 km/s, which corresponds roughly to the velocity spac- 
ing adopted for the searching grid used to correct AV„/ M (see 
above), which implies that this velocity grid was well adapted 
for our purpose. 



250 



200 



.150 



01 

> 



100 - 



50 - 




50 100 150 200 
AV™ S ' [km/s] 



250 



Fig. 2. Comparison between the AV mo dei values obtained in 
Monte-Carlo simulations of 100 GIRAFFE data-cubes using 
a 0.8 and 1 arcsec seeing (see text). The black line is a linear 
fit fixing the intercept to zero. The residual dispersion is ~ 4.3 
km/s. 



To find the best model, we finally looked for the (r t ,Vf[ at ) 
pair minimizing the difference between AV m ,d e i and AV b s - We 
checked that, in all simulations, such a criterion gives a unique 
solution. Results are listed in TableQ] and representative exam- 
ples of this kinematical fitting for three RD galaxies are shown 
in Fig. [3] 



3.2. Testing the method of correction 

To test the correction method, we performed Monte-Carlo sim- 
ulations of 100 GIRAFFE data-cubes. The parameter space in- 
vestigated encompassed the inclination, half-light radius, PA, 
rotation velocity, and RC gradient, of values randomly cho- 
sen from the typical values observed for the GIRAFFE sam- 
ple. Half of the Monte-Carlo simulations were generated using 
an arctan RC shape model, and the other half using an expo- 
nential term combined with a second-order polynomial term. 
This polynomial term allowed us to create a central velocity 
"bump" as observed in some early-type galaxies (see Fig. [TJ. 
We note, for late-type RCs, that the Vfiat did not have to be 
located within the IFU, which allowed us to test how the cor- 
rection procedure operated in this case. For early-type galax- 
ies, Vfiai should be located within the IFU, as V mnx is gener - 



. g., Courteau 
(IFlores et alJ 



1997) 



2006). 



ally located close to 2.2Rd ~ \315Ri ia if (e.: 
which almost always falls within the IFU 
We note, however, that in some cases, because of the dynamical 
influence of the bulge, the RC of some early -type galaxies can 
show an extended veloci ty peak up to large distances (see, e.g., 
M3 1 . [Carignan et al . 2006). The limited spatial coverage of the 
GIRAFFE IFU would probably provide us with an overesti- 
mation of the rotation velocity. In the RD sample, all galaxies 
besides one have B/T < 0.2, which means that such an effect 
cannot affect significantly the RD sample. 

These simulated RCs were used to simulate GIRAFFE ob- 
servations following the method outlined in Sect. 3.2. These 
simulated data-cubes were in turn used as inputs to test the 
method of correction on AV i, s . In Fig. [4] we compare the input 
rotation velocity Vfiat, with the asymptotic velocity Vo of the 
best arctan model obtained using the method detailed above. 
We find a good linear correlation with Vo = [0.65 ± 6.22] + 
[1.00 + 0.03] V// a ,, consistent with Vo = Vfi at , and a residual 
dispersion ~ 17 km/s. In this figure, we have distinguished be- 
tween simulations where Vfi at is sampled by the IFU (black 
circles) from those where this is not the case (black squares). 
If we compare the corrected velocity Vo with the last point of 
the RC sampled by the IFU V e „d, both sets of simulations fall 
on the same region. In such a plot (not shown here), we find a 
similar result, with V = [-1.23 + 4.47] + [1.03 + 0.02] V md , 
also consistent with Vo = Vend- 

We checked that r, had little influence on the AV b s correc- 
tion by comparing the corrected rotation velocity, derived by 
allowing r, to vary between the three values defined in Sect. 
3.1, with that obtained by fixing r, to the middle value of this 
grid, i.e. the one that makes most of the RC gradient fell inside 
one GIRAFFE pixel as observed in most cases. We found that 
both sets of values correlated well within ~ 6 km/s (1 - cr), i.e., 
as expected, r, has a negligible impact on the derived correc- 
tion. This is due mainly to the large seeing size (~0.8 arcsec) 
compared to velocity gradients associated with typical r, values 
(see Fig.HJ, which significantly dilutes the variations of the RC 
gradient, once projected onto the IFU. 

It is important to check whether or not these simulations 
are representative of real galaxies. To d o this, we selected a 



few galaxies from the GHASP survey (lAmram et alJ 12002). 
which provides us with high spatial-resolution data-cubes for 
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J033219. 68-275023. 6 



z=0. 55955 Inc. =58 RD 



J033212.39-274353.6 



3 



z=0.42130 Inc.=90 ED 



J033234. 04-275009.7 



z=0. 70162 Inc.=59 RD 




Fig. 3. Examples of kinematical fitting of three z ~ 0.6 RD galaxies. From left to right: HST/ACS F775W image with the 
GIRAFFE IFU superimposed (from Paper I), observed VF (shown with a 5x5 interpolation see Paper I), best modeled VF 
(shown with a 5x5 interpolation), residual map between the observed and modeled VFs. Relatively large differences are found 
only close to the minor axis, w here departure from pure circular motion is artificially exaggerated by projection effects (see a 
discussion of this effect in, e.g.. lChemin et al ] |2006l) . 



a morphologically-complete sample of local galaxies using 
Fabry-Perot interferom etry. The analysis of the GHASP sam- 
ple is still underway (lEpinat et al.ll2.008D . and a full compari- 
son between local 3D data and distant 3D GIRAFFE data will 
be addressed in a future work. Here, we check whether or not 
the simulations are representative of local galaxies, especially 
in the range of velocity and size spanned by distant RDs. We 
restrict our choice to galaxies for which an arctan shape pro- 
vides a relatively good representation of the RC, because we 
probe the accuracy of the correction applied to AV b s , and not 
to the measurement itself, which is an independent issue. Given 
the limited number of galaxies meeting these criteria, we chose 
to restrict our choice to 7 such galaxies with rotation veloci- 
ties ranging from 190 to 220 km/s, where are typical measure- 
ments for most distant RDs (see next Sect. 4). This choice al- 
lows us to sample roughly the most relevant velocity range for 
distant RDs, rather than testing a few isolated points spread 
over the full velocity range of the distant sample. We sim- 
ulated GIRAFFE observations by degrading the Fabry-Perot 
high-resolution data-cube to the resolution of GIRAFFE obser- 
vation (0.8 arcsec seeing) and then to the GIRAFFE IFU spatial 
sampling. From these data-cubes, we extracted a VF and AV bs 
as for real GIRAFFE data, and corrected them. These simula- 
tion results are shown as blue stars datapoints in Fig. [4] which 
agree with the Monte-Carlo simulations that fall in the same 
velocity range, reproducing the dispersion of the correlation. 
This confirms that we can confidently rely on the correction 
applied to AV b s - 

In Fig. [5] we plot all correction factors a = Vfi a ,/AV b s 
obtained for the Monte-Carlo simulations. Both simula- 
tions of redshifted local galaxies and real GIRAFFE rotat- 



ing disk galaxies fall in the same region of the plot, ap- 
part from three galaxies (J033212.39-274353.6, J033230.78- 
275455.0J033245. 11-274724.0), that show relatively large a. 
One is a compact galaxy (J033245.1 1-274724.0, see Paper 
II), while another one is seen almost edge-on (J033212.39- 
274353.6), which might explain the relatively high a. The last 
galaxy (J033230. 78-275455.0) is a special case, because only 
half of the galaxy is detected in emission (this galaxy was 
shifted into the RD+ class in Sect. 2. 1). We checked that the re- 
sults presented in this paper are not significantly affected when 
these three special objects are removed from the sample. In the 
RD subsample, we found a median a ~ 1 .26, consistent with 
what is found in the Monte-Carlo simulations, with a mean 
a= 1.25 ±0.12. 

3.3. Random uncertainty budget of the rotation 
velocity derivation 

The uncertainty budget can be decomposed into uncertainties 
related to AV bs> an d those related to the correction applied to 
AVobs to obtain V/i at . 

The main source of uncertainty that can affect AV bs is that 
associated with a finite spectroscopic SNR in the measurement 
of Vj™ and Vj^ , which can be quantified using Monte-Carlo 
simulations (see Fig. [6]). For each galaxy, we used the SNR 
maps derived in Paper I to estimate the mean SNR uncertainty 
on Vj££ and VjSF, i.e., on AV b s - We found a median (mean) 
uncertainty due to a finite SNR of ~9 km/s (8km/s). 

We used Fig. [4] to quantify the uncertainty associated with 
the correction made to AV (> bs- We found 17 km/s (1 — a resid- 
ual dispersion), with no noticeable systematic uncertainties. 
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Fig. 4. Comparison between the V/uu used as inputs to the 
Monte-Carlo simulations of 100 GIRAFFE data-cubes, with 
the Vo values obtained after using the method of correction de- 
scribed in Sect. 3.2. Black dots represent RCs generated using 
an arctan model (i.e., late-type like RCs), while open circles 
represent RCs typical of early-type galaxies. Black squares rep- 
resent simulations where the rotation velocity is not sampled 
by the IFU (see text). The black line is a linear fit, and the 
dash lines represent the 1 - cr residual dispersion ~\1 km/s. 
Blue stars represent real observations of local galaxies artifi- 
cially redshifted to z ~ 0.6 (see text). 



We note that this uncertainty takes into account cases where 
the RC is not fully sampled by the IFU. As we have already 
pointed out, one galaxy in the sample deserves more consider- 
ation. Only part of the galaxy J033230.78-275455.0, which was 
classified as RD+ (see Sect. 2.1), was detected in emission (see 
Paper I). AV b s does not encompass a symmetric range along 
the RC, which biases the rotation-velocity input of the model 
to a lower value. Considering a rough RC mo del, composed of 
a first , linear, rising part up to R max = 2.2Ro dPersic & Salucci 
199 lb . and a second flat part, we estimate a ~30km/s correction 
to "symmetrize" AV /, S . This correction is convolved directly 
by the other uncertainty factors on AV b s - Following Sect. 3.1 
and 3.2, we also took into account uncertainties due to seeing 
and r, variations (~5 and 6 km/s, respectively), which produced 
a total uncertainty for the correction on AV h s of ~19 km/s. 

The total uncertainty on V// fl ,xsin (z) (i.e., the radial rotation 
velocity uncorrected for inclination) was derived by adding in 
quadrature all the previous terms. We finally propagate the un- 
certainty associated with the inclination, which is ±5 degrees, 
according to Paper I. We finally found a median (mean) total 
uncertainty on V fiat of ~37 km/s (44km/s). The 1 - cr disper- 
sion of the total uncertainties around the mean is ~30 km/s. 



1.5 - 




R half [ 



Fig. 5. Correction factors a used to correct AV b s for RDs 
(blue dots; open blue dots represent the two RD+ galaxies). 
Small black dots correspond to Monte-Carlo simulations of 100 
GIRAFFE data-cubes in the same range of half-light radius, 
PA, inclination, rotation velocity, and RC gradient (see text). 
The horizontal dash line represents to the mean correcting fac- 
tor of 1.25 derived from Monte-Carlo simulations, while the 
1 - cr dispersion around the mean is shown in dash lines. 
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Fig. 6. Monte-Carlo simulations of the velocity measurement 
accuracy (see Paper I for details). The red line shows the cor- 
responding bias (almost zero), while the black dashed lines 
shows the 1 - cr error: ~12 km/s between SNR=3-5, ~5 km/s 
between SNR=5-10, and < 2 for SNR>10. 



4. The K-band Tully-Fisher relation at z ~ 0.6 

4.1. Results 

We compare both the distant and local K-band TFR in Fig. 
[7] Holding the slope to the local value, cr res increases from 
RDs (0.31 mag), to PRs (0.80 mag), and CKs (2.08 mag). 
Restricting the distant TFR to dynamically well-relaxed RDs, 
the local and distant relations have comparable ow we there 



fore confirm quantitatively the previous findings of Flor es et al 
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(2006) that all of the enlarged dispersion of the distant TFR 
comes from non-relaxed galaxies. If we allow the slope to vary 
during the fit, we find a residual dispersion of 0.31 mag. This 
indicates that assuming no evolution in slope appears to be a 
reasonable assumption. 

The zero point of the TFR of distant RD galaxies is found 
to be -5.88 ± 0.09 (1 — cr bootstrapped error, slope fixed), i.e., 
0.66 mag fainter that the local zero point. Allowing the slope 
to vary during the fit, we find a zero point of -5.07 + 2.37 and 
a slope of -7.24 + 1.04 in agreement (within the correspond- 
ing uncertainties) with those derived by fixing the slope to the 
local value. This strongly suggests no evolution in slope of the 
K-band TFR. Given the limited number of RDs in the distant 
sample (16), we adopt the zero-point value derived by keep- 
ing the slope to the local value. We note that if we consider 
the RD+ class (18 galaxies, see Sect. 2), we find similar results 
(see Tabled. 

We find in the distribution of Mk residuals, a skewness 
and a kurtosis of -0.22 and -0.99, respectively for RDs: this 
is roughly consistent with the Gaussian residuals, within the 
1 - cr expected thresholds^. This indicates that residuals are not 
biased significantly on any side of the relation. A Welcher's t- 
test then provides a probability <scl% that the local and distant 
relations have the same zero point. We note that during the fit- 
ting procedure, we weighted all the rotation velocities by their 
associated uncertainties. Therefore, the result of the Welcher's 
t-test means that within random uncertainties, the zero-point 
difference of 0.66 mag between the distant and the local rela- 
tions is statistically significant. 




log(V flat ) [km/s] 

Fig. 7. Evolution of the K-band TFR (AB magnitudes). The 

completude limit Mk 20.14 (corresponding to Mj = -20.3, 

see Sect. 2) is indicated by an horizontal dash line. Blue dots 
represent RDs (the two RD+ galaxies are represented with 
open blue dots), green squares PRs, and red triangles CK galax- 
ies. The black line is the local TFR, while the blue dash-line 
represent a linear fit to the z ~ 0.6 TFR (see text). 



4.2. Systematic uncertainty budget 

We investigate which systematic effects could bias the evolu- 
tion of TFR zero point detected between z ~ 0.6 and z = 0. 



4.2.1. Photometric systematic uncertainty 

As stated in Sect. 2, we do not expect any noticeable 
systematic-effect on measurements due to photometric mea- 
surement or calibration. However, one possible exception is 
the extinction correction made on both local and distant galax- 
ies. Even if we follow exactly the same procedure for the two 
samples, one might wonder whether or not the low surface- 
brightness outer-regions in distant galaxies could be far less 
well detected compared to local galaxies, which would create a 
bias in the derivation of their axis ratios and extinction correc- 
tions. However, we found in Paper II that the HST/ACS imag- 
ing used to derive the axis ratio in distant galaxies allows us to 
reach the optical radius (~ 3.2 x Rj), which limits such a bias. 
We note that extinction corrections are relatively small in the K- 
band: in the distant sample, we find a mean/median correction 
of 0. 1 1 ±0.05 mag ( 1 - cr dispersion), while in the local sample, 
the mean (median) correction is found to be 0.15(0.14)±0.05, 
which are consistent within the respective uncertainties. 



4.2.2. Effects related to the comparison to the local 
relation 

One important aspect about the evolution of the TFR is to con- 
trol systematic effects between the distant and local relations. 
These considerations ha ve led us to adopt a s a local refer- 
ence the TFR derived by Hammer et al. (2007) for a complete 



subsample of the SDSS dPizagno et alj|2007l) (see Sect. 2.2). 
However, there are several differences between the local and 
distant samples that could lead to possible biases when com- 
paring the corresponding TFRs. 

First, kinematical data are obtained using different tech- 
niques (3D spectroscopy vs. long-slit spectroscopy). We miti- 
gated this effect by restricting the lo cal sample to galaxie s hav- 
ing a well-defined RC (see details in Hammer et al.l2007 ). This 
also mitigates the assumption that both local and distant galax- 
ies have RC shapes than can be well described by an arctan 
(see Sect. 2.2 and 3.2). Second, a different emission line is 
used for deriving the rotation velocity ([Oil] and H a ). However, 
both emission lines are f ound to provide similar estimates of 
the rotation velocity (e.g.. lSaintongd l2007). Third, in the local 
sample, Vgo (defined as the velocity at a radius that encom- 
passes 80% of the light) is used as a proxy for V/; a ,, while 
we directly used the asymptotic Vo from the arctan RC in the 
distant sample . Vso i s a reliable proxy for Vf i at, as shown in 
Hammer et al.1 (120071) (see also IPizagno et al.ll2007l) . Finally, 
one possible caveat about the use of 2MASS K-band magni- 
tudes was recently pointed out by Noordermeer & Verheijenl 



2 For a distribution of n points, these thresholds are V(6/'i) and 
V(21/«) respectively (n=16 is the present case). 
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Table 2. Fits to the local and distant K-band TFRs, using M K {AB) = a + b x log (V/; n ,). 



K-band TFR 


Slope a 


Zero point b 


Comments 


Local relation 


-6.88±0.57 


-6.54+1.33 


SDSS subsample 


Distant RDs 


-6.88 


-5.88+0.09 


Using local slope 


Distant RDs 


-7.24±1.04 


-5.07±2.37 


Slope free 


Distant RD/RD+ 


-6.88 


-5.92±0.10 


Using local slope 


Distant RD/RD+ 


-6.47±1.30 


-6.85±2.98 


Slope free 



( 120071) . i.e. that 2MASS underestimates the K-band luminos- 
ity in relatively low surface brightness galaxies. We note that 
this bias would go in the opposite trend compared to the evo- 
lution seen in the TFRs, i.e., to shift the local relations to- 
wards brighter magnitudes or larger stellar masses. However, 
we have also mitigated this effect by restricting the local rela- 
tion to galaxies having Log(Vfi ar ) > 2.2, i.e. to the most mas- 
sive galaxies (see Sect. 2.2). 

4.2.3. Kinematic systematic uncertainty 

We restricted the analysis of the TFR evolution to well-relaxed 
RDs. In distant TFR studies, it is often assumed that the (gas) 
rotation velocity can be simply derived as the half-amplitude of 
the VF divided by sin (/) (modulo spatial-resolution effects): an 
important assumption behind this simple prescription is that the 
azimuthal velocity component of the 3D velocity vector within 
the galaxy (i.e., the rotation) dominates its projection along 
the line-of-sight. This assumption should always be checked 
a posteriori, by showing that residuals between a rotating-disk 
model and the observed VF are small compared to the ampli- 
tude of rotation motions. In the case of GIRAFFE RDs, there 
are few doubts that this assumption is correct, because their 
modeled VFs match the observed ones (see Paper I). 

Convers ely, CK galaxies ar e clearly out of dynamical equi- 
librium (see lFlores et alj|2006l and Paper I). As a consequence, 
one does no longer know which component of the 3D velocity 
vector dominates its projection along the line-of-sight. If these 
galaxies are associated with mergers, as we will discuss below, 
deriving the rotation velocity in this way is probably uncertain, 
if not meaningless. This is the reason why in this study, we 
restricted the analysis of the TFR to well-relaxed RDs. 

In Sect. 3, we tested extensively the method of correction 
used to correct AV t, s for spatial-resolution effects. All galaxies 
were assumed to be RDs, reg ardless of their kinematical class 
(see iPuech et a l. 2006, 2007). However, this provides uniform 
a values that are independent of kinematical class: a misclas- 
sification of a galaxy has no impact on the way a is derived. 
This approach helps to compare with long-slit spectroscopy re- 
sults, where all galaxi es are implicitly assumed to be RDs (e.g., 
Conselice et al.ll2005l) . As explained above, the fit to the TFR 
is restricted to RDs only, and therefore this cannot affect the 
evolution of the zero point of the TFR. 

One might wonder whether or not the slight bias found in 
Fig. |4] between early- and late-type galaxies could influence 
significantly the results. Once translated into log ( V// fl/ ) (using 
Eq. [TJ, the offset of the distant TFR is found to be -0. 1 dex 
between z ~ 0.6 and z = 0, which is much larger than that 



found in Fig.|4]for early-type galaxies (~0.025 dex). Late-type 
galaxies have the opposite trend, which provides a -0.014 dex 
offset in Fig. [4] It is impossible to explain all of the shift in the 
distant TFR zero point in terms of such a morphological bias. 
Looking at the TFR residuals against B/T for RDs analyzed in 
Paper II, we find no special trend, which excludes the presence 
of a bias. 

Another source of systematic uncertainty could be the lim- 
ited spatial coverage of the GIRAFFE IFU. Fig.|4](see the black 
squares) shows that when the plateau of the RC is not sampled 
by the IFU, the recovered rotation velocity is underestimated 
on average by ~0.03 dex. We note that this effect generally 
leads to the underestimation of V fiat, which would increase the 
shift of the zero point between the local and the distant TFRs. 
However , most of the RDs are well spatially covered by the 



IFU (see iFlores et al.l l2006 and Paper I), and this could affect 
a few galaxies in our sample. It is thus unlikely that such an 
effect could affect significantly the results in a systematic way. 

Finally, the most important possible systematic uncertainty 
likely comes from our limited knowledge of the seeing dur- 
ing observations. Individual variations from galaxy to galaxy 
leads to relatively small random uncertainties (see Sect. 3.2). 
However, we assumed a uniform value of 0.8 arcsec during 
the rotation-velocity correction process. According to Fig. [2] 
changing the seeing in the simulations from 0.8 to 1.0 arcsec 
implies a systematic effect of -0.02 dex on the rotation velocity 
correction. 

4.2.4. Total systematic uncertainty 

The only possible systematic effect we can identify so far is 
that associated with the kinematics, due to the correction for 
the rotation velocity. Using Eq.Q] this translates into a possible 
systematic uncertainty of ±0. 14 mag in the 0.66 mag evolution 
of the K-band TFR zero point between z ~ 0.6 and z — 0, which 
represents 20% of the shift. 

5. Discussion 



5.1. Comparison wit mFlores et all (2000) 

In IFlores et all d2006l) (see also IPuech et"aD I2007I) . we used 
hydro-dynamical simulations of an Sbc galaxy to infer a mean 
correction factor a - 1.2 + 0.05 on V max for galaxies of typical 
diameters between 2 and 3 arcsec: this value is less than that 
found here using Monte-Carlo simulations. At first sight, it is 
surprising to find a larger mean a for V/; fl , than for V mn x (since 
V fiat < Vmax). The ex planation is likely linked to the fact that 



in lFlores et al. d2006l) . we did not consider all possible ranges 
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of size, PA, inclination, rotation velocity and RC gradient to 
derive this mean correction, because we were limited by the 
hydro-dynamical models of spiral galaxies available to us at 
that time. Another differen ce is the exponenti al dependence of 
a with galaxy size found in lFlores et al . ( 2006), which is notre - 
produced here. The reason for this is that in lFlores et a l. ( 2006), 
we used the same simulation of an Sbc galaxy to simulate dis- 
tant galaxies of different sizes, which induced an intrinsic cor- 
relation between the RC gradient and the galaxy size. Our tests 
show that if we introduce this correlation in the Monte-Carlo 
simulations, we then recover the exponential variation of a with 

Rhalf- 

The bri ghtening found in th e K-band TFR was not initially 
detected by Flores et al. ( 20061). In Fig.[8l we sh ow the K-band 



-26 



TFR obtained following the lFlores et alj (2006) method to cor 



rect AV bs, i.e., appl ying a cons t ant co rrection factor of 1 .2. The 
black line shows the lVerheiienl d200ll) K-band local TFR using 
Vfi a t as a kinematical tracer for th e rotation velocity, i .e., the re- 
lation used as a local reference by lFlores et al.l (12006). Both the 
local and distant relation are then in good agreement: if we fix 
the slope of the distant relation to the local one, we find a shift 
~ 0. 1 mag between the zero point of the distant and th e local 
relations. Therefore, we retrieve the F lores et al. (2006) results 
that no significant evolution in zero point can be detected. The 

reason for such a difference is twofold. 

d2007l) . 



First, as detailed in Hammer et al 



crepa ncy between the K-band TFR derived in the Verheiien 
(2001) sampl e, compared to those derived in t wo other lo- 
cal sa mples dCourteaull 19971 Piza gno et al.ll2007l) . We show in 
Hammer et al. ( 20071) that this discrepancy is due to the larger 
fractio n of faint and slow rotating galaxies in the Verhe iien 
(2001) sample compared to the two other ones. Such galax- 
ies have a larger gas fraction, which results in a different 



slope at the low-mass e n d of the TFR dMcGa ugh 2005). 
Noordermeer & Verheijenl (12007) derived a new K-band lo- 
cal TFR: they found a slope similar to th e one infe r red by 
Hammer et al. (2007) after r estricting th e Verheiier] (2001) 
TFR to the high-mass end. In Flores et alL 2006b. the rela tion 



used as a local reference was however the I Verheiien ( 200 lb 



lation. If we refit the distant RD subsample usi ng the rotation 
velocities derived following iFlores et al. ( 2006b . but using the 
slope of the SDSS local relation (i.e., the one used in Fig. [7] 
see the blue lines in Fig. [HJ, we find a shift ~0.4 mag between 
the d istant and the local zero points. Compared to Flores et all 
(2006), we use a different local sample, which provides a better 
control of systematic effects that can occur in the comparison 
with the distant sample (see Sect. 4.2.2): this rigorous approach 
allows us to explain ~ 60% of the TFR zero-point shift, previ- 
ously hidden by spurious effects. 

Second, we have significantly improved the method for 
correcting the rotat ion velocity since the preliminary work of 



Flores et alj (|2006). This is illustrated in Fig. |9] where we 



compa re the rotation velocities derived following Flo res et al 



(2006) and t hose d erived in this study: this figure reveals that 
Flores et al. (2006) underestimated the rotation velocity in RDs 
by 1 1 % on average, which corresponds to a ~0.05 dex shift in 
log (Vfiat), or ~0.3 mag once converted into a Mk shift using 
Eq. Q] Therefore, we attribute the remaining ~40% of shift in 




2 2.5 
l°g(V max ) [km/s] 



Fig. 8. K-band TFR derived following the lFlores et"al] d2006b 
methodol ogy, for the RD subsample. The black line is the local 
relation oflVerheiienl ( 200lb . i.e., the one used as a reference by 
Flores et al.l (12006). Note that K-band magnitudes in the dis- 
tant sample have been converted into the Vega system using 
MxiVega) = Mk(AB ) - 1.85. For sim plicity, we assume that 
the K' magnitudes of I Verheijenl (1200 lb are roughly equivalent 
to those derived in the distant sample using the ISAAC Ks fil- 
ter. The blue line is the SDSS local TFR used as a reference in 
this study (converted into the Vega system). The blue dash line 
is a linear fit to the distant RDs, which has a zero point 0.4 mag 
lower than the local one, fixing the slope to the local value. 
Open symbols represent galaxies from the new CDFS sam- 



ple, w hile full symbols represent galaxies used in Flo res et al 
(200i). 



the TFR zero point, previously undetected, to the improvement 
in the rotation-velocity derivation. 

Finally, Fig.[8]provides us with a useful cross-check of the 
new method used in this study to correct the rotation velocity. 
In Fig. [8] one RD galaxy (J033212.39-274353.6) appears to be 
shifted to lower velocities compared to other distant rotating 
disks. This galaxy was identified as having a particularly high 
a value in Sect. 3.3, while its data are well described by the 
TFR in Fig. [7] This suggests that the new method of correc- 
tion appears to be better suited for deriving rotation velocities: 
comparing Fig.[8]with Fig. [7] we note that the dispersion of the 
RD subsample is significantly reduced using the new method 
of correction. 

5.2. Origin of the scatter in the distant relation 

Using the spatially-resolved kinematics provided by 3D spec- 
troscopy, we confirm earlier results that the dispersion of the 
distant relations can be ex plained by the pres ence of dynami- 
cally non-relaxed galaxies dFlores et al. 2006b . These galaxies 
are found to be out of dynamical equilibrium, with random mo- 
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Fig. 9. Comparison betwee n the rotation vel ocities VV06 ob- 
tained using the method of Flore s et al. I d2006l) (i.e., a constant 
correction factor of 1.2) vs. rotation velocities Vfuu derived us- 
ing the new method used in this paper, for the subsample of 
RDs. Wo6 underestimates Vfyt by ~\\% on average. Open 
symbols represent galaxies from the n ew CDFS sample, w hile 
full symbols represent galaxies used in lFlores et al. U2006I) . 



tions instead of ordered motions, especially for the most com- 
pact galaxies ( Puech et alj|2006l ; IWeiner et al.l l2006). Ongoing 
comparisons between hydrodynamical simulations of galaxy 
mergers with GIRAFFE dat a reveal that they could be asso- 
ciated with major mergers dPuech et alJl2007l : Puech et al., in 
prep.). Such events provide a natural link between CK galax- 
ies and morphologically -peculiar galaxies (see Paper II), which 
causes an increase i n the dispersion of the TFR at high redshift 
dKassin et alJl2007l:lAfkinson et al.|[2007b . 

On the other hand, it is still unclear whether or not more 
quiescent processes, such as secular gas ac cretion through in 



ternal processes (bars), or cold-gas flows (Dekel & Birnboim 
20061: Kere s et al 1 l2005h . could explain such anomalous kine- 
matics on a relatively large spatial scales, especially for CK 
galaxies (see Paper I). More work is required to compare dy- 
namical predictions of such events with observations. 



5.3. Comparison with the evolution of the K-band 
luminosity density 

Finding a 0.66 mag brightening of the K-band TFR between 
Z ~ 0.6 and z = appears to be quite surprising, given the 
opposite trend seen in the evolution of the K-band luminos- 
ity density in "blue" galaxies over the same redshift range. 
One may reasonably assume that "blue" star-forming galax- 
ies and emission-line galaxies belong to th e same population 

d2007l) find that the K- 



( Ham mer et alj 11997ft : lArnouts et al. 



band luminosity density of blue star-forming galaxies fades by 
0.5-0.6 mag from z ~ 0.6 to z = 0, which is exactly the op- 



posite trend. What appears to be a clear contradiction simply 
reflects two different methodologies. 

As we have detailed above, the TFR allows us to compare 
two physically-connected populations of galaxies, i.e., the lo- 
cal spirals and the distant rotating disks: the evolution of this 
relation directly reflects the evolution of the K-band luminos- 
ity in rotating disks of similar total mass between z ~ 0.6 and 
z — 0, assuming tha t the rotation velocit y can be used as a proxy 
for the total mass dHammer et alj 120071) . On the other hand, 
the color selection used to select "blue" galaxies at z = and 
at z ~ 0.6 can produce heterogeneous galaxy populations, as 
noticeable when comparing their luminosity densities. Optical 
colors are well known to be seriously affected by instantaneous 
star formation and extinction. The evolution of the K-band lu- 
minosity density in "blue" galaxies reflects the number evolu- 
tion of "blue" galaxies, w hich are much more numerous at high 
redshift (e.g.. lEllis|[l997l) . 

To check that there is no contradiction between the evo- 
lution of the TFR zero point and the evolution of the K-band 
luminosity density, we derived the K-band luminosity den- 
sity in the complete GIRAFFE sample, i.e., including all dy- 
namical classes. As stated in Sect. 2.1, the GIRAFFE sam- 
ple is representative of z ~ 0.6 emission-line galaxies with 
Mj(AB) < -20.3 (i.e., with M steUar > 1.5 x 10 10 M o , see Paper 

I) . We therefore expect a K-band luminosity density in this 
sample that represents the luminosity density of "blue" galax- 
ies at these redshifts. We estimated the K-band luminosity in 
the GIRAFFE sample to be log(p^) = 8.74 L .Mpc 3 , which 
agree well with the K-ban d luminosity density of "blue" star- 
forming galaxies found by Arnouts et al.l (2007). This compar- 
ison confirms that the GIRAFFE sample is representative of 
z ~ 0.6 emission-line galaxies: the K-band luminosity density 
in the GIRAFFE sample is consistent with the results inferred 
from studies using far larger samples. 

5.4. Did distant rotating disks double their stellar-mass 
over the last 6 Gyr? 

Once restricted to well-relaxed RDs, we find, between z ~ 0.6 
and z = 0, a shift in the TFR zero point of 0.66±0.14 mag in 
Mk, or -0. 1 +0.02 dex in log (V// a r). We consider the interpreta- 
tion of this shift in terms of galaxy evolution. 

Could there be any so-called "progenitor bias" between dis- 
tant RDs and local relaxed spirals, which could imply that the 
latter would not be the descendants of the former? Local RDs 
are found to be twice as numerous as distant ones (see Paper 

II) , which implies that some were not in a relaxed dynami- 
cal state at z ~ 0.6 (being then PR or CK). However, local 
intermediate-mass spirals have only a low probability to have 
undergon e a major merger since z ~ 0.6 (15-30%, as dis- 
cussed in Puech et alj|2007l) : this means that most distant RDs 
must dynamically evolve smoothly towards local relaxed spi- 
rals, which implies that distant RDs are the progenitors of a 
majority of local spirals. Therefore, it makes sense to interpret 
the evolution of the RD-restricted TFR as an evolution in lu- 
minosity, rotation velocity, or a combination of the two. We 
reiterate that we assume no evolution in slope (see Sect. 4.1): a 
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larger sample would be required to tackle directly the possible 
evolution in slope of the TFR and determine the details of the 
evolutionary path between distant and local RDs. 

We discuss the possibility that this shift could correspond 
to a pure luminosity brightening of 0.66+0. 14 mag in RDs over 
the past 6 Gyr. Such a luminosity-brightening in RDs would 
correspond to a growth in stellar-mass, which can be estimated 
in the following way. Between z ~ 0.6 and z — 0, the evo- 



lar mass for med since z=l occu rred in star-forming galaxies 



range 

0.16 dex dProry et alj2004tlArnouts et alJ2007l) . depending on 
the selection criteria. Using \og{M steUar ) = log(M ste u ar /L K ) + 
log (Lk), one finally finds a stellar-mass evolution of 0.39-0.42 
dex. A more exhaustive derivation, using the stellar-mass TFR, 
leads to a similar conclusion, with an evolution in zero point 
of 0.36° ( , 1 06 dex between z ~ 0.6 and z = (see Appendix 
A). On the other hand, the evolution of the cosmic stellar-mass 
densi ty over the same redsh i ft range is found t o be ~0.1-0.16 
dex (iDickinson et al J 120031, iDrorvet alJ l2004t lArnouts et al 



20071) . iHammer et al.1 (2005) found that most of the present- 
day stellar-mass forme d since z = 1 occur in intermediate- 



mass galaxies (see also lBell et al. 



2005), which include 60% of 



emission-line galaxies at z ~ 0.6 (Ham mer et al.ll 1997b . Hence, 
if we assume that, on average, quiescent galaxies do not evolve 
in stellar-mass, one can estimate the growth in stellar-mass be- 
tween z ~ 0.6 and z = in intermediate-mass, emission-line 
galaxies to be ~0. 32-0. 38 dex. Therefore, this interpretation 
of the TFR zero-point shift as a pure luminosity evolution is 
roughly in agreement with the evolution of the cosmic stellar- 
mass density over the same redshift range. This interpretation 
is supported by the fact that distant RDs are found to be LIRGs 
or at least star-bursting, i.e., they are actively forming stars at 
very high rates (see Paper II). 

Interpreted in this way, the evolution of the zero point 
of the K-band TFR reflects the growth in stellar-mass of the 
most active population over the past 6 Gyr, i.e. star-forming 
intermediate-mass galaxies, by a factor ~2.5-2.6. This com- 
pares well with the gaseous O/H phase abundance of z ~ 
0.6 emission-line galaxies, which is, on average, half that 
of present-day sp irals for a similar range of stellar mass 
( Liang et alj20 06). Such an evolution in stellar mass would im- 
ply that RDs converted an important amount of gas into stars 
over the past 6 Gyr. In other words, the main evolutionary path 
for RDs during the past 6 Gyr would be conversion of gas into 
star through gas supply, which is further supported by their 
relatively low gas-disk V/cr values JPuech et alj|2007l) and the 
inside-out build-up of their stellar disks (see a detailed discus- 
sion in Paper II). 

The opposite interpretation is that the shift in the zero point 
of the K-band TFR corresponds to a pure rotation-velocity evo- 
lution of -0. 1 +0.02 dex between z ~ 0.6 and z = 0. Even in this 
case, the observed evolution in M ste ii ar /Lfc over this redshift 
range still implies, on average, a substantial growth in stel- 
lar mass of 0.13-0.16 dex in intermediate-mass galaxies (see 
above). As a consequence, the stellar mass growth in quies- 
cent intermediate-mass galaxies should be approximately simi- 
lar because the stellar mass growth in the entire population (i.e., 
quiescent or not) is observed to be similar in this redshift range. 
This is clearly not what is observed, because ~80% of the stel- 



(LIRGs, see IHammer et al. 2005). We can therefore rule out a 
pure evolution along the velocity axis. 

The last possibility is a combined evolution along both 
axes, i.e. a simultaneous brightening in luminosity with a de- 
crease in rotation velocity V rot . According to the Virial theo- 
rem, Vj ot scales as the ratio of the total mass enclosed within 
the optical radius R opr , over R opt (~ 3.2Rd in both distant and 
local galaxies, see Paper II). A decrease in V rot between z ~ 0.6 
and z = would then imply a decrease in this ratio, which in 
turn would imply that R op , (or R c /) increases faster than the to- 
tal mass over the same redshift range. Between z ~ 0.6 and 
z = 0, Rd does not seem to evol ve strongly, at lea st in the RD 
subsample (see Paper II; see also lPuech et alF0 07). Therefore, 
only a moderate increase in the total mass enclosed within R opr 
could occur over this redshift range. Such a scenario would 
agree with observations: as stated above, local intermediate- 
mass spirals have a low probabilit y to have undergon e a major 
merger since z ~ 0.6 (15-30%, see lPuech et al ■1 l2007l and refer- 
ences therein). Therefore, most distant RDs should be already 
mostly assembled at z ~ 0.6, and should not undergo strong 
evolution in terms of total-mass from z ~ 0.6 to z = 0. 

To explore the real evolution along both axes, we com- 
pare our results with t he model of spherical gas accretion of 



Birn boim et alJ (l2007h . lBirnboim et alJ (120071) present a model 
of accretion for a star-bursting galaxy at ~0.7 and describe 
its subsequent evolution in terms of mass, down to z = 0. 
Although it is not clear whether or not such a model could 
be representative of all properties of distant RDs, we assume 
that it can be used to constrain the average ma ss evolution in 



a typical RD halo. In that case, the results of iBirnboim et al 



(2007) would suggest that the baryonic mass in the disk re 
mains constant between z ~ 0.6 and z = 0, while the Virial 
baryo nic mass roughly do ubles (see their Fig. 2). On the other 
hand. lConroy et al.1 (120071) found that the Virial-to-stellar mass 
ratio in intermediate-mass galaxies is roughly constant between 
z = 1 and z = 0. Assuming that the su bsample of RDs fol lows 
the same trend, one can combine the Conrov et al. ( 2007|) ob 



servational results with the model of Birnboim et al 



(E007l) . 



and find that RDs would roughly double their stellar-mass be- 
tween z ~ 0.6 and z = 0. Accounting for the evolution in mass- 
to-light ratio over this redshift range, this would translate into 
an evolution of 0.35-0.43 mag in luminosity, to be compared 
with the 0.66 mag found in the evolution of the TFR zero point. 
Once translated into log (V// af ), this would allow a -0.04 dex 
evolution along the velocity axis, between z ~ 0.6 and z — 0. 

In conclusion, we find that the most likely interpretation of 
the evolution of the TFR zero point is the one in which this shift 
reflects mostly a luminosity evolution of RDs. We estimated an 
upper limit to the contribution of an evolution along the veloc- 
ity axis to be at most one half of the total shift of zero point. 
Such a brightening of distant RDs over the past 6 Gyr would in- 
dicate a doubling of their stellar-mass, which is independently 
supported by their other dynamical and morphological proper- 
ties. 

Finally, this could suggest that the baryonic (stars plus gas) 
TFR would not evolve with redshift: if one accounts for the 
(average) two times larger gas fraction of z ~ 0.6 galaxies 
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compared to z — dLiang et al. 120061) . one finds that dis- 
tant galaxies roughly fall back onto the local smTFR (see 
Appendix A). Studies of the local baryonic TFR have shown 
that galaxies having V/i at <100 km/s systematicall y fall be- 
low th e TFR defined by more massive galaxies (e.g.,|McGaugh 
20051) . However, once the gas fraction is accounted for, all 
galaxies follow the same baryonic TFR. Interpreted that way, 
our results suggest an evolution of this threshold toward higher 
masses (velocities) at higher redshifts, of at least Vflat ~300 
km/s at z ~ 0.6. This supports the idea that the baryonic TFR 
would be much m ore "fundamental" that the stellar-mass TFR 
dMcGaughl2005l) . 



6. Conclusion 

We have studied the evolution of the K-band TFR, using a 
representative sample of 65 emission line, intermediate-mass 
galaxies atz ~ 0.6, unaffected by field-to-field variations within 
Poisson statistics, and observed using 3D spectroscopy. We 
have presented and tested a new method allowing us to safely 
recover V/iat of distant galaxies. We have also investigated pos- 
sible sources of systematic effects, which can be particularly 
important when studying the evolution of the TFR. We have 
paid particular attention to the analysis of both the local and 
distant samples in as similar as possible ways. We draw the 
following conclusions: 

(1) The larger dispersion of the distant TFR is caused by 
galaxies with anomalous kinematics, ranging from per- 
turbed rotators to very dynamically-disturbed galaxies. We 
find a positive and strong correlation between the complex- 
ity of the kinematics and the scatter contributed to the TF. 
Once restricted to relaxed rotating disks, the scatter of the 
TFR appears to not have evolved, which might suggest no 
evolution in slope; 

(2) We detect for the first time a significant evolution of the 
K-band TFR zero point, which we attribute to an average 
brightening of z ~ 0.6 galaxies by 0.66±0.14 mag. We at- 
tribute this evolution to the fact that selected distant galax- 
ies are either starbursts or LIRGs. The distant emission-line 
rotating disks represent roughly one fourth to one fifth of 
present-day spiral progenitors, and one half of the whole 
population of z ~ 0.6 intermediate-mass rotating disks (see 
Paper II). Therefore, a significant part of spiral progenitors 
are doubling their stellar masses during the last ~6 Gyr, in 
good agreement with for mer studies on the e volution of the 
mass-metallicity relation dLiang et al.l l2006). which would 
suggest no evolution of the baryonic TFR with redshift; 

(3) Current studies of the evolution of the TFR, even using 
spatially-resolved kinematics, are limited by important sys- 
tematic uncertainties, which can be attributed to the limited 
spatial resolution of the kinematics, and to the derivation of 
the stellar-mass. These systematic uncertainties represent 
20% of the evolution in zero point of the K-band TFR. 

Further progress in the study of the evolution of the TFR 
will benefit from a forthcoming larger sample of RDs. This will 
allow us to reduce the random uncertainties, allow the slope to 



vary during the fit, and start studying TF residuals. Although 
we have limited the local sample to galaxies having a well- 
defined rotation curve, a decisive answer about the evolution of 
the TFR will require both local and distant representative sam- 
ples processed in the same way, from the observational strat- 
egy (3D spectroscopy) to the data analysis itself, as we have 
done in this paper. A similar analysis of a representative lo- 
cal sample observed using 3D spectroscopy is required to con- 
firm the above results. This point will be addressed in a future 
study, using the local kinematic GHASP survey, whose analysis 
is currently ongoing dEpinat et al. 2008). Finally, an important 
limitation remains linked to systematic uncertainties, which are 
mostly due to a lack of spatial resolution. The only way to over- 
come this limitation will be to develop a new generation of op- 
tical integral-field spectrographs with higher spatial-resolution, 
i.e., assisted by adaptive optics. In this respect, the future ex- 
tremely large telescopes will allow us to make a decisive le ap 
forward in our understanding of the TFR dPuech et al . 2008). 



Appendix A: The evolution of the stellar-mass TFR 

A.1. The local stellar-mass TFR 

Stellar masses M s , e u „ r were estimated from M ste u ar /LK ratios 
using the method of Bell et alJ ([2003). We used a solar lumi- 



nosity in the Ks-band of 3.28 (Vega) and assumed a "diet" 
Salpeter IMF dBell et al.ll20"ol . This method takes advantage 
of the tight correlation found between rest-frame optical colors 
and M s , e n ar /Lic ratios, assuming a universal IMF. These cor- 
relations are found to be relatively insensitive to the details 
of ga l axy SFH, dust c ontent, and metallicity dBell & de J ong 
200 It iBell et alj|2003l) . which implies that they are invaluable 
for deriving stellar mass without being too sensitive to the de- 
tails of the stellar population synthesis models. We note that 
following this method, M ste n ar /LK ratios are corrected for the 
amount of light d ue to red-g i ant st ars using g-r colors. 

According to IBell et alJ d2003l) . using this method, the to- 
tal random uncertainty on log(M s , e u ar /M Q ) at z ~ should 
be lower than 0. 1 dex, and the systematic uncertainties due to 
galaxy ages, dust, or bursts of star-formation can reach 0. 1 dex. 
Finally, the influence of TP- AGB stars in the derivation of stel- 
lar masses cou ld result in an oyerestimation of the stellar mass 
by -0.14 dex (Marasto n et al.ll2006l: IPozzetti et al] E007) in a 



systematic way. 

Using this method, we converted the local K-band TFR in 
the subsample of the SDSS (see Sect. 2.2) into stellar masses 
and found: 

log(M stellar /M ) = 4.46±0.53 + (2.8+0.23)xlog(V /to ),(A.l) 
with 0-^5=0.15 dex. 

A.2. The distant stellar-mass TFR 

To estimate stellar masses in the distant sample (see Table [TJ, 
we used the same method as in the local sample, i.e., the K- 
band luminosity and B-V rest-frame colors. One important is- 
sue is whether or not the same correlations between color and 
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stars. However, they found that this systematic is constant with 
redshift, at least up to z ~ 1.2. In other words, at a given 
IMF, relative predictions between two different redshifts of SPS 
models are much more robust in terms of M/L predictions. 

We therefore focus in Fig. IA.2I on the evolution of 
log (M ste u ar / ' L K ) between z ~ 0.6 and z - 0: we show the his- 
tograms of log {M stellar I Lk) derived using the lBell et a l. (2003) 
method, both in the distant and local samples. Both histograms 
were centered using the median log {M s teiiar I Lk) of the local 
sample: we find that log (M ste u ar / Lg) evolves from z ~ 0.6 to 
z — by ~0.06 dex. In this figure, we show other determina- 
tions of the evolu tion of log (M „,,» „ , /L y ) from the literature. 
Drorv et alj d2004l) . using lMarastonl i 1998b SPS models, found 
an evolution of 0. 13 dex in galaxies with stellar masses between 
4xlO lo M and 10 n M o , i.e., in a range of stellar-mass similar 
to GIRA FFE galaxies. This is similar to the 0. 15 dex evolution 



2 2.5 
log(V flat ) [km/s] 



Fig. A.l. Evolution of the stellar-mass TFR in the RD subsam- 
ple (the two RD+ galaxies are represented with open blue dots). 
The black line is the local smTFR, while the blue dash-line rep- 
resent a linear fit to the z ~ 0.6 smTFR 



found by Arnouts et alJ d2007l) in a flux-limited samp le over the 
same redshift ran ge, usin g Bruzua l & CharloJ d2003l) SPS mod- 



Arnouts et al J d200~ ) find an evolution of 0. 1 6 



M ste ii ar /LK ratios found at z=0 by Bell et al. (2003) can be di 



els. Furthermore, 
dex once restricted to a sample of blue star-forming galaxies. 
Because such blue galaxies and emis sion-line galaxies pr oba- 
bly belong to the same populations dHammer et alJ[l997l) . we 
conclude that we might be underestimating the evolution of 
log (M ste u ar j Lk) by up to ~0. 16-0.06=0. 1 dex in a systematic 
way. 

What is the origin of this systematic effect? As 
log {M stella r I Lk) depends mostly on color and not on mass, 
et al J (12003b used SPS models to predict their 



Bell . 



' aver- 



rectly applied at higher redshift. iBorch et al.N 2006l) showed a S e " correlations given a reasonable range of SPS parameters 



that this appears possible at least up to z~0.6 (see their Fig. 
4), with an associated random uncertainty of ~0.3 dex, and an 
average systematic uncertainty of up to -0.2 dex (i.e., towards 
an overestimation of the stellar-mass at high z). 

We show the derived distant smTFR in Fig. I A. 1 1 The 
smTFR shows the same trend with the kinematic classification 
as the K-band TFR. Restricting ourselves to RDs, and holding 
the slope constant, the distant and local relations have simi- 
lar dispersions (cr rcj =0.12 dex). If we allow the slope to vary 
during the fit, we find a residual dispersion of 0.12 dex, sug- 
gesting again that assuming no evolution in slope is accept- 
able. Maintaining the slope at its local value, the distant smTFR 
zero point is found to be 4. 10±0.16(random) ^(systematic, 
see above), i.e., 0.36 dex smaller than the local zero point (see 
Table lA.lT i. A Welcher's t-test gives a probability <k1% that the 
two relations have the same zero point. 



A. 3. Stellar-mass systematic uncertainties 

The most important systematic uncertainty is that associated 
with the mass-to-light ratios [M/L] predicted by stellar popu- 
lation synthesis [SPS] models. Absolute values of M/L depend 
mostly on the SPS model and IMF used, and more particularly 
on the prescript i ons for the TP-AGB s tellar-evolution phase 



(Maraston 2005; Maraston et al. 2006). Pozzetti et al. (2007) 



com pared M/L pred ictions between iBruzual & Chari ot (2003) 
and Maraston (2005) SPS models and found a systematic dif- 
ference of -0.14 dex due to different prescriptions for TP-AGB 



(e.g., m etallicity, star formation histories). However. iBell et al 
(2003) did not explicitly fit the age but assumed instead a ref- 
erence age of 12 Gyr at z = 0. This introduces a systematic 
differenc e with other r esults where age is fitted explicitly. For 
instance, iDrorv et alj ([2004) find an average age of 8.1 Gyr 
at z = 0. If one compares the median log {M s teiiar I Lk) found 
in the local SPSS sample used in th is study, with the aver- 
age value found by lDrory et al. (2004) in a similar mass range, 
one finds a difference of -0.08 dex (which accounts for the dif- 
ferent IMF used). Similarly, we find a difference of -0.16 dex 
between the median log {M ste u ar / Lk) fo und in the G IRAFFE 
sample, and the average value found by IDrorv et aD d2004l) at 
z = 0.6. These values are consistent with the systematic error 
bars quoted above. They suggest that we underestimate the evo- 
lution of log (M s t e n ar / ' L K ) between z — and z ~ 0.6 by -0.08 
dex, which is consistent with Fig. IA.2l 

In summary, these comparisons show that we are probably 
underestimating the evolution of log (M stellar I 'Lk) between z 
and z ~ 0.6 by -0.08 dex, and maybe up to -0.1 dex due to th e 
use of the simplified prescriptions for M/L o flBell et al.1 (120031) . 
This directly translates into a systematic effect of the evolution 
of smTFR zero point by up to +0. 1 dex. 

A.4. Total systematic uncertainty 



We summarize all possible sources of systematic uncertainty 
identified so far in Table IA.2I We expressed all in terms of 
their influence on the shift in the zero point between the dis- 
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Table A.l. Fits to the local and distant smTFRs, using log (M ste u ar I M Q ) — a + b x log (Vfi at ). 



K-band TFR 


Slope a 


Zero point b 


Comments 


Local relation 


2.8±0.23 


4.46±0.53 


SDSS subsample 


Distant RDs 


2.8 


4.10+0.16 


Using local slope 


Distant RDs 


3.28+1.10 


3.01±2.47 


Slope free 


Distant RD/RD+ 


2.8 


4.13+0.15 


Using local slope 


Distant RD/RD+ 


2.61+1.33 


4.56±3.03 


Slope free 
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Fig.A.2. Histograms of log (M ste iiar I Lk) found in the local 



and di stant samples of galaxies using the method of iBell et al 



(2003). Both histograms have been re-centered using the me- 
dian found in the local samples, which allows us to directly in- 
fer the evolution of log (M ste u ar /L K ) between z ~ 0.6 and z — 0, 
i.e., ~0.6 25 dex. Also shown , the evolution of log (M ste u ar /L K ) 
found by Drorv et all d2004l) in a sam ple of intermediate-m ass 
galaxies (black long-dashed line), and [ Arnouts et al.l (12007b in 
blue star-forming galaxies (blue mixed-line) or independently 
of the color (black mixed-line). 



tant and the local smTFRs. Taking into account all these un- 
certainties, we find a shift in the smTFR zero point of 0.36_qL 
dex between z ~ 0.6 and z — 0. Systematic uncertainties rep- 
resent up to ~60% (in dex) of the evolution detected in zero 
point, i.e., much larger than systematic uncertainties affecting 
the evolution of the K-band TFR zero point. The derivation of 
stellar mass is highly model-dependent, which implies a large 
additional uncertainty. However, we note that these effects tend 
to overestimate the stellar mass at high redshift, i.e., to min- 
imize the evolution of the smTFR. Using log {M ste u ar I Lk) = 
log (Medlar) - log (L K ), one finds that log (M ste „ ar IL K ) evolves 
by -0.1 dex (= 0.36 - 0.66/2.5) between z ~ 0.6 and z = 
in the subsample of RDs. Therefore, the evolution in the zero 
point of the smTFR is quite robust, in the sense that we might 



have minimized the evolution of M ste tt ar /Lfc through the choice 
of method used to derive the stellar mass (see Sect. A. 3). 
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